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Abstract 


Theoretical  and  experimental  studies  were  carried  out  for  the  transport  of  cations  in  the 
channel  between  a polymer  coating  and  a metal  substrate  from  a defect  in  the  absence  of  an 
electrical  applied  potential.  The  model  consists  of  two  stages:  an  initial  period  during  which 
ions  diffuse  in  the  coating/metal  interfacial  "channel"  and  adsorb  on  the  coating  surface  and  a 
propagation  period  during  which  ions  also  diffuse  into  the  coating.  The  mathematical  models 
were  solved  to  predict  the  cation  concentration  and  flux  under  the  coating  and  the  relative  rate 
of  diffusion  between  initial  and  propagation  periods.  Model  parameter  values  were  derived  from 
the  results  of  an  experiment  conducted  in  a specially  designed  diffusion  cell.  The  experiment 
measured  the  depletion  of  Na"^  ions  in  a cylindrical,  central  reservoir,  which  was  placed  within 
the  perimeter  of  a defect  scribed  through  the  coating  of  an  epoxy-coated  steel  panel.  Model 
predictions  of  concentration  versus  time  agreed  well  with  experimental  results,  which  showed 
that  most  of  Na"^  ions  were  removed  by  lateral  diffusion  from  the  reservoir  during  the  initial 
period.  Further,  the  transport  during  the  initial  period  was  much  faster  than  that  during  the 
propagation  period.  The  results  also  indicated  that  during  the  propagation  period,  the  rate 
limiting  step  was  the  lateral  diffusion  along  the  coating/metal  interface  rather  than  diffusion 
through  the  coating. 
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1.  Introduction 


Corrosion  of  metals  beneath  organic  coatings  is  a problem  of  major  importance  to  coating 
technologists.  In  order  for  corrosion  to  occur  there  must  be  sources  of  oxygen,  water,  and  an 
electrolyte  to  transfer  charge.  When  the  corrosion  occurs  beneath  the  coating,  these  chemicals 
can  be  transferred  by  diffusion  or  an  applied  or  induced  potential  to  the  corrosion  site.  Any  of 
a number  of  processes  may  control  the  overall  corrosion  rate,  but  diffusion  of  chemical  species 
is  always  a likely  candidate  to  limit  transport  because  it  is  inherently  a slow  process.  Typical 
diffusion  rates  through  coatings  are  orders  of  magnitude  less  than  their  unhindered  values. 

Migration  of  cations  can  occur  by  diffusion  between  anodic  and  cathodic  corrosion  sites 
at  the  interface  between  the  coating  and  the  metal.  These  sites  may  be  initially  close  to  one 
another,  but  with  time  they  can  separate  and  become  localized.  Cathodic  blisters  can  form  and 
eventual  coating  failure  may  result.  The  transport  of  cations  such  as  sodium  provides  a 
mechanism  to  neutralize  the  hydroxide  ions  (1): 

O2  + 2H2O  + 4e-  + 4Na-"  = 4(Na+OH-) 

In  the  presence  of  an  applied  or  induced  potential,  diffusion  is  facilitated.  It  is  possible 
to  produce  concentrations  within  cathodic  blisters  which  are  many  times  greater  than  their  values 
in  the  external  solution.  This  has  been  observed  experimentally  and  predicted  theoretically  (2). 
The  model  of  Pommersheim  et  al.  (2)  involves  two  periods:  an  initiation  period  and  a 
propagation  period.  During  the  initiation  period  ions  diffuse  through  the  channel  between  the 
coating  and  the  substrate.  An  advancing  front  where  cations  are  instantaneously  removed  by 
sorption  is  predicted  to  move  between  the  anode  and  the  cathode.  During  the  propagation 
period,  cations  break  through  into  the  cathodic  blister  where  they  accumulate.  During  both 
periods,  the  applied  cathodic  potential  facilitates  the  diffusion  process.  A critical  value  of  the 
potential  was  predicted  below  which  blisters  would  not  grow  and  above  which  they  would. 

In  this  paper  the  transport  of  cations  (sodium  ions)  in  the  channel  between  an  organic 
coating  and  a metal  substrate  is  studied  in  the  absence  of  an  applied  electric  field.  Both 
diffusion  along  the  coating/ metal  interface  and  diffusion  outward  through  the  coating  (in  the  axial 
direction)  are  considered  as  well  as  adsorption  of  sodium  ions  on  the  underside  of  the  coating. 

2.  Experimental 

2.1.  Specimen  Preparation 

Figure  1 presents  a schematic  of  the  specimen  configuration  for  the  measurement  of  the 
diffusion  of  cations  at  the  organic  coating/steel  interface.  It  consists  of  an  organic  coated 
substrate,  an  inner  glass  cylinder  containing  a solution  of  sodium  chloride  and  a large,  outer 
cylinder  holding  distilled  water.  The  substrate  was  a SAE  1010,  cold-rolled  steel,  and  the 
organic  coating  was  a commercial,  two-component,  room-temperature  cure,  pigmented  epoxy 


1 


resin.  The  total  pigment  content  (by  mass)  in  the  coating  was  42.5  percent.  The  substrate  had 
a matte  finish  and  surface  roughness  ranging  from  0.9  to  1.3  ^m.  The  substrate,  having  a 
dimension  of  152x102x0. 8mm,  was  cleaned  by  repeatedly  and  alternately  rinsing  with  acetone 
and  methanol  (reagent  grades),  followed  by  dry  air.  The  coating  was  applied  on  the  substrate 
immediately  after  drying,  using  the  draw-down  technique.  Excess  coating  was  placed  on  one 
end  of  the  substrate  panel  then  pulled  down  along  the  length  of  panel  using  a glass  rod. 
Masking  tape  strips,  5mm  wide,  placed  along  the  length  of  the  panel  edges  were  used  to  control 
the  coating  thickness.  The  coating  was  allowed  to  cure  at  room  conditions  (24 °C  and  45% 
relative  humidity)  for  two  weeks.  The  thickness  of  the  coating  was  approximately  150  /xm  as 
determined  by  an  eddy  current  gage. 

After  curing,  a disk  of  the  coating  13mm  in  diameter  was  bored  to  the  steel  substrate  and 
removed  from  the  center  of  the  coated  panel.  A glass  tube  (inner  cylinder),  having  an  outside 
diameter  of  12.5  mm,  an  inner  diameter  of  12  mm,  and  a height  of  50  mm,  was  placed  where 
the  coating  disk  was  removed  (Figure  1).  Epoxy  adhesive  was  used  to  seal  the  space  between 
the  glass  tube  and  the  coating.  An  0.05-mm  thick  gasket  made  of  tetrafluoroethylene  polymer, 
having  an  outside  diameter  of  13mm  and  an  inside  diameter  of  10  mm,  was  situated  at  the 
bottom  of  the  glass  tube.  This  procedure  ensured  that  the  epoxy  adhesive  did  not  come  into 
contact  with  the  steel  substrate,  which  could  have  blocked  the  lateral  diffusion  route  of  cations. 
Using  a room-temi>erature-cure  silicone  adhesive,  an  open-ended. 


Figure  1.  Specimen  configuration  for  measuring  the  diffusion  of  cations  beneath  an  organic 
coating  on  a substrate. 
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PMMA,  outer  cylinder,  with  a height  of  50  mm  and  an  inner  diameter  of  50  mm,  was  attached 
to  the  coated  substrate,  symmetrically  around  the  inner  cylinder.  The  water-filled,  outer  cylinder 
was  used  so  that  the  coated  specimen  was  completely  soaked  with  water. 

2.2.  Cation  Measurement 

The  diffusion  of  cation  in  the  "channel"  between  the  coating  and  the  metal  substrate  was 
measured  by  following  the  decrease  of  Na"^  ion  concentration  in  the  inner  cylinder.  The 
concentration  of  Na"^  ions  in  the  inner  cylinder  was  measured  using  a commercial,  1-mm  tip, 
Na'^-miniselective  electrode.  The  electrode  was  connected  to  a high-impedance  (>  10^‘*  ohm) 
ion  analyzer.  The  ion  analyzer  was  programmed  to  record  Na"^  ion  concentrations  at  desired 
intervals.  The  ion-selective  electrodes  respond  to  activity  rather  than  concentration,  where 
activity  is  related  to  the  concentration  by  the  activity  coefficient.  However,  because  the 
measurement  system  was  calibrated  in  terms  of  concentration  we  reported  our  results  in  terms 
of  concentration.  The  Na'^-selective  electrode  was  calibrated  prior  to  the  actual  measurements 
and  its  reliability  was  ensured  by  checking  at  the  conclusion  of  the  experiment,  A preliminary 
experiment  indicated  that  in  order  to  maintain  the  same  liquid  volume  in  the  inner  cylinder 
throughout  the  measurement,  the  coated  substrate  must  be  saturated  with  water.  For  that  reason, 
double-distilled  water  was  placed  in  the  outer  cylinder  24  hours  prior  to  filling  the  inner  cylinder 
with  the  cation-containing  solution,  a 0.4  moles/L  solution  of  NaCl  in  double-distilled  water. 
Thus,  it  is  expected  that  there  was  a water  layer  at  the  coating/substrate  interface  within  the  area 
covered  by  the  outer  cylinder  prior  to  filling  the  inner  cylinder  with  the  cation  solution  (3). 

The  top  of  the  inner  cylinder  was  sealed  with  molten  wax  and  cellophane  film  was  used 
to  cover  the  top  of  the  outer  cylinder  (Figure  1).  This  procedure  was  used  to  minimize  the 
evaporation  of  water  from  both  cylinders.  The  duration  of  the  experiment  was  174  hours. 
There  was  little  change  in  the  liquid  levels  in  both  the  inner  and  outer  cylinders  during  this 
period.  There  was  some  corrosion  on  the  bare  steel  area  inside  the  inner  cylinder.  Further,  at 
the  end  of  the  experiment,  the  coating  was  easily  detached  from  the  substrate,  indicating  that  the 
adhesion  was  completely  lost  after  the  experiment.  However,  the  metal  underneath  the  coating 
was  still  bright  and  showed  no  evidence  of  corrosion.  The  data  collected  together  with  the 
specimen  characteristics  are  presented  in  the  Appendix  I. 

3.  Conceptual  Model 

The  sodium  ions  diffuse  from  the  inner  chamber  (Figure  1)  into  the  "channel"  which  lies 
between  the  coating  and  the  metal  substrate.  There  are  no  concentration  gradients  of  sodium 
ion  within  the  cylinders.  The  transport  process  is  presumed  to  occur  in  two  stages  or  periods. 
During  the  first  period,  ions  diffuse  through  the  channel  to  a point  (rj)  where  they  are  rapidly 
removed  by  adsorption  onto  the  bottom  surface  of  the  coating.  This  process  defines  a moving 
interface  (located  at  r;),  which  advances  under  the  coating.  The  initiation  period  ends  at  time 
0;  when  the  interface  reaches  the  far  end  of  the  coating  (r=R),  During  the  second  or  last 
period,  which  begins  at  time  0;,  diffusion  under  the  coating  still  occurs,  but  this  time  the  sink 
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for  the  ions  is  diffusion  of  ions  into  the  coating  from  the  channel.  Thus,  in  this  period  diffusion 
occurs  in  two  directions,  along  the  coating/metal  interface  (radial  direction)  and  normal  to  the 
coating  (axial  direction).  This  process  is  assumed  to  occur  by  a first  order  process,  so  that 
where  the  concentration  of  ions  in  the  channel  is  higher,  more  diffusion  into  the  film  will  occur. 
During  both  periods  ions  are  being  removed  from  the  solution  in  the  inner  cylinder,  but  during 
the  first  period  ions  are  removed  more  rapidly  since  diffusion  distances  are  smaller.  Thus,  the 
second  stage  should  occur  over  time  periods  longer  than  the  first. 

Since  higher  concentrations  of  sodium  are  present  near  the  entrance  of  the  channel 
(at  r=a),  more  transport  into  the  coating  will  occur  there.  The  amount  of  adsorption  and 
diffusion  into  the  coating  will  depend  on  the  thicknesses  of  the  channel  and  coating,  the  cell 
dimensions  and  the  diffusivity  of  sodium  ions  through  the  channel  and  coating.  These  factors 
will  also  determine  whether  transport  is  controlled  by  diffusion  along  the  coating/metal  interface 
or  through  the  coating. 

4.  Mathematical  Model 

4.1.  Model  I (Initial  Stage) 

During  the  initial  stage  the,  A-cations  (considered  to  be  sodium  ions  here)  are  diffusing 
under  the  influence  of  a concentration  gradient.  The  concentration  of  sodium  ions  is  0^(0  at 
r=a  and  zero  at  r=ri  where  they  are  removed  by  rapid  adsorption.  This  establishes  a linear 
concentration  gradient. 

The  molar  flow  of  A,  w^  is  given  by: 

0 

dC,  c 

w . = -2  7c  r 6 D — - = -2  IT  6 D^— 

^ ^ dr 

r 

where  a is  the  radius  of  the  inner  cylinder,  r is  any  radius,  6 is  the  height  of  the  channel,  and 
De  is  the  effective  diffusivity  of  sodium  ions  in  the  channel. 

At  steady  state  the  diffusive  flow  must  be  equal  to  the  rate  at  which  sodium  ion  is  adsorbed 
or  bound  to  the  surface: 


2 n 6 DC„ 

e a 

]n^ 

a 


(1) 


dr. 

' — - 2 n r. 
^ dt 
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where  Ps  is  the  surface  molar  density  (mol/cm^)  of  sodium  ions  in  the  adsorbed  layer. 
Equating  Eqs.(l)  and  (2)  gives: 


-P. 


dt 


5 D,  CJ,f) 


(3) 


The  concentration  of  sodium  ions  Ca(t)  within  the  inner  chamber  can  be  found  from  an  ion 
balance  at  r=a: 


= ^A\r=a^  ~^e' 


ac, 

dr 


^ 2 n a 6 

= a 


(4) 


where  V is  the  volume  of  solution  in  the  inner  cylinder. 

The  sodium  ion  concentration  can  be  found  as  a function  of  interface  position  by  accounting 
for  the  total  amount  of  adsorbed  sodium  ion: 
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where  C^o  is  the  initial  sodium  ion  concentration  in  the  inner  cylinder  and  C is  the  dimensionless 
concentration. 

Substituting  Eq.(5)  into  Eq.(3): 
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where  y is  the  dimensionless  distance  q/a,  and 
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a,  the  reciprocal  of  the  time  constant  for  the  initial  period,  is  given  by: 

^ ^ ^ (8) 
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Eq.(5)  and  (6)  can  be  used  together  to  predict  the  concentration  of  sodium  ion  within  the 
inner  cylinder  as  a function  of  time.  The  initial  stage  ends  when  t = 0;  at  y = R/a.  For  fixed 
values  of  the  parameters  a and  y,  Eq.(6)  can  be  used  to  predict  0;.  As  indicated  by  Eq.(7),  y 
is  a dimensionless  parameter  which  is  a measure  of  the  amount  of  adsorption  which  can  occur. 

4.2.  Model  n (Final  Stage) 

The  second  or  final  stage  begins  when  the  first  stage  ends  (at  t=0i).  During  this  stage, 
diffusion  occurs  both  along  the  coating/metal  interface  and  into  the  coating.  This  model  was 
developed  by  first  making  a sodium  ion  balance  on  a thin  shell  in  the  channel  located  at  radius 
r and  having  width  Ar: 


2Tzr&  NX-2n  & 


2 Ji  r Ar  D/C^  - 
h 


(9) 


where  5 is  the  height  of  the  channel,  is  the  radial  flux  of  sodium  ions,  D,.  is  the  diffusivity 
of  sodium  ions  through  the  coating,  h is  the  coating  thickness  and  C^w  is  the  concentration  of 
sodium  ions  in  the  outer  cylinder. 

In  the  limit  as  Ar  0 we  obtain: 
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So  that 
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(11) 


In  the  above  formulation,  several  assumptions  are  implicit.  The  flux  of  sodium  ion  into 
the  coating  is  taken  equal  to  Dc(CA-CAw)/h-  This  formulation  presumes  that  a linear 
concentration  gradient  of  cation  across  the  coating  is  rapidly  established  and  maintained,  a 
situation  which  is  made  more  likely  because  the  coating  is  thin.  In  addition,  the  diffusivities 
and  Dg  are  taken  constant.  This  is  a reasonable  assumption  considering  the  uniformity  of  the 
coating  and  the  fact  that  the  free  stream  or  bulk  diffusivity  is  relatively  independent  of 
concentration  for  sodium  chloride  solutions  (4).  Also,  since  the  amount  of  water  in  the  outer 
cylinder  is  over  ten  times  that  of  the  solution  in  the  inner  cylinder  (Appendix  I),  will  remain 
low  as  compared  to  C^. 

Eq.(ll)  is  subject  to  the  boundary  conditions: 


C^=C^(0  at  r=a 


and 


=0  at  r=R 


dr 


The  last  condition  arises  because  there  should  be  little  or  no  Na"^  ion  migrating  beyond  the 
end  of  the  channel  (at  r=R)  (underneath  the  wall  of  the  outside  cylinder).  The  relatively  small 
area  here  (2x6R)  insures  that  only  a small  number  of  ions  will  be  transported  there  as  compared 
to  those  that  pass  into  the  much  larger  area  beneath  the  coating. 

Eq.(lO)  and  the  boundary  conditions  can  be  made  dimensionless  by  the  substitutions: 


e = 


(12) 


and 


(13) 
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leading  to: 


^ ^ dV 


h b D, 


Q = Q 


with  dimensionless  conditions: 


and 


e = 1 at  i = - 
R 


at  I = \ 


(14) 


Alternatively,  Eq.(14)  can  be  written: 


^2  ^ _ ^2  ^2  0 ^ Q 


(15) 


where 


2 R^ 
= — 

Dhb 

e 


(16) 


The  dimensionless  parameter  m is  a measure  of  the  relative  importance  of  diffusion  along 
the  coating/metal  interface  to  diffusion  through  the  coating.  If  m is  large,  diffusion  along  the 
interface  will  control  transport,  while  if  m is  small,  diffusion  through  the  coating  will  control. 
Because  the  geometric  ratio  R^/h5  will,  in  general,  be  quite  large  for  a typical  coating  system, 
if  the  diffusion  coefficients  are  of  comparable  magnitude,  diffusion  along  the  interface  is  more 
likely  to  control  transport. 
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The  solution  of  Eq.(15)  is  given  in  terms  of  modified  Bessel  functions: 

e = Cj  (m  O + q (m  0 


07) 


Applying  the  boundary  conditions,  the  constants  Cj  and  C2  can  be  found,  so  that: 

(m)  (m  0 + /j  (m)  (m  O 
y — 

i:,  (m)  /,  (m|)  + /,  (m)  A;  (m-|) 

/v  /v 


Eq.(18)  governs  the  dimensionless  concentration  profile  within  the  channel  during  the  final  stage. 

As  shown  by  Eq.(4),  a sodium  ion  balance  on  the  inner  cylinder  equates  the  depletion  of 
sodium  ions  to  the  amount  which  diffuse  out  through  the  channel  entrance  (at  r=a).  The 
gradient  ACJAr  at  r=a  is  found  by  differentiating  Eq.(18).  Performing  this  operation  and 
simplifying: 


2 D C 

1 c a 


(19) 


where  the  parameter  rj  is  given  by: 


2 af  hP^h 


(20) 


with 


^ /j  (m)  (m  a/R)  - (m)  /j  (m  a/R) 

“ AT,  (m)  (m  a/R)  + /,  (m)  K„  (m  a/R) 

Eq.(19)  is  integrated  to  give  the  sodium  ion  concentration  as  a function  of  time: 
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(22) 


where  C’^ois  the  concentration  of  sodium  ion  at  the  start  of  the  final  stage  (t=0i). 
a\  the  reciprocal  of  the  system  time  constant  for  model  n,  is  given  by: 


a'  = 

Vh 


2 n a b r\ 


(23) 


Eq.(22)  predicts  that  the  concentration  will  fall  exponentially  for  times  greater  than  0;. 
Together  with  Eqs,(5)  and  (6),  which  are  valid  for  times  less  than  0;,  it  governs  the  way  that 
the  concentration  of  sodium  ions  within  the  inner  cylinder  vary  with  time.  Model  constants  can 
be  determined  by  fitting  experimental  data  to  these  equations. 

5.  Results  and  Discussion 

The  data  given  in  the  Appendix  I for  the  sodium  ion  at  different  times  within  the  inner 
cylinder  can  be  used  to  evaluate  the  best  values  of  the  model  constants  a and  y (for  model  I)  and 
a’  (for  model  II).  For  the  range  of  experimental  data  reported,  y must  lie  between  0 and 


0.0765. 


The  initiation  time  0;  depends  on  both  a and  y since  from  Eq.(6): 


R 


a 


(24) 


a,  a’  and  7 were  found  by  minimizing  the  overall  sum  of  squares  of  the  deviations  between  the 
observed  and  predicted  concentrations  for  both  models.  A numerical  simulation  was  used  to 
evaluate  the  integrals  in  Eq.(6)  and  Eq.(24).  The  computer  programs  for  these  operations  are 
given  in  Appendix  11. 

The  best  values  of  the  model  constants  for  this  method  were  found  to  be 
O'  = 9.5  X 10'^  h'^  Q!’=  5.0  X 10'^  h'^  and  7 = 0.065.  Using  these  values,  0;  was  evaluated 
as  49  hours  by  using  Eq.(24).  Figure  2 is  a plot  which  compares  the  original  data  to  that 
predicted  by  the  models.  A good  fit  is  obtained  for  model  11,  while  for  model  I,  for  times  less 
than  0i,  the  fit  is  not  as  good.  This  disparity  may  arise  because  of  unsteady  state  effects 
associated  with  the  initial  penetration  of  sodium  ions  into  the  channel.  Both  models  assume  that 
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diffusion  occurs  at  quasi-steady  state,  so  that  concentrations  within  the  channel  have  enough  time 
to  continually  readjust  to  the  changing  concentrations  within  the  inner  cylinder.  However,  there 
may  be  a small  initial  period  during  which  quasi-steady  state  is  not  yet  established.  The  epoxy 
coating  will  take  a short  time  to  become  acclimated  to  the  first  sodium  ions  which  enter  the 
channel.  During  this  period,  rapid  adsorption  may  occur.  Diffusion  into  the  coating  can  also 
be  slowed  by  the  accumulation  of  sodium  ions  in  the  distilled  water  in  the  outer  cylinder  so  that 
the  assumption  that  Caw=0  will  not  be  valid  at  longer  times.  Unfortunately,  no  sodium  ion 
concentration  data  were  collected  in  the  outer  cylinder  to  verify  this  assumption.  However,  even 
if  no  adsorption  occurred,  the  highest  value  of  Caw  could  only  reach  688  ppm  under  the 
experimental  conditions  employed.  With  adsorption  removing  the  bulk  of  sodium,  the  highest 
value  of  Caw  would  fall  to  107  ppm. 


Figure  2.  Comparison  of  predicted  to  experimental  values  of  sodium  ion  concentration  at  the 
defect  as  a function  of  time  (C^  vs  t). 
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The  ratio  of  the  constants  a/a’  is  a measure  of  the  time  constant  Tq/ti,  where  Tq  is  the  time 
constant  for  model  II  and  Tj  is  the  time  constant  for  model  I.  For  the  experimental  data  obtained 
(Appendix  I),  this  ratio  is  1.9,  indicating  that  the  initial  diffusion  and  adsorption  of  ions  on  the 
channel  occurs  on  a time  scale  faster  than  their  subsequent  diffusion  into  the  coating.  The  value 
of  7[(R/a^-l]  is  the  amount  of  sodium  ion  that  can  adsorb  in  the  channel  relative  to  the  amount 
that  is  potentially  available  for  adsorption.  For  the  experiment  conducted  this  is  0.85,  and  from 
Appendix  I the  matching  concentration  occurs  at  about  42  hours.  This  indicates  that  the  coating 
adsorbs  all  but  15%  of  the  sodium  ion  available  within  that  time.  The  remainder  is  then 
available  to  diffuse  into  the  coating.  Note  that  the  models  predict  that  this  will  occur  in  49 
hours. 

From  the  best  values  of  a and  y,  the  diffusion  coefficient-channel  depth  product 
Dg5  is  1.71  X 10’^  cmVs.  The  ease  with  which  the  epoxy  coating  could  be  detached  from  the 
steel  substrate  at  the  conclusion  of  the  experiment  suggests  that  5 may  be  large.  For  an  average 
channel  width  of  6=0.1  mm,  this  would  give  an  effective  lateral  diffusivity  of  0^=  1.71  x 10'^ 
cm'Vs,  which  is  close  to  the  free  stream  value  for  sodium  chloride,  1.5  x 10'^  cm'Vs  (4)  and 
three  times  higher  than  that  of  Na"^  ions  through  a sodium  exchange  membrane  (5). 

Figure  3 presents  a plot  of  the  parameter  r]  versus  the  modulus  m.  a/R  is  shown  as  a 
parameter,  rj  is  a quantitative  measure  of  the  penetration  of  ions.  It  is  the  ratio  of  the  molar 
flow  (moles/s)  of  ions  beneath  the  coating  to  the  maximum  flow,  which  would  occur  if  the  whole 
underside  of  the  coating  could  be  exposed  to  solution  having  the  external  concentration  C^o.  It 
is  a dimensionless  measure  of  how  effective  the  coating  is  in  removing  cations.  values  lie 
between  0 and  1.  Small  values  of  m yield  large  rj  values  (~  1),  while  at  large  values  of  m,  rj 
becomes  small.  With  the  maximum  flow,  there  are  no  concentration  gradients  beneath  the 
coating  and  the  modulus  m is  small.  As  shown  by  the  figure,  rj  is  then  near  unity.  In  general, 
when  m < 0.5  diffusion  through  the  coating  will  control  transport,  while  if  m > 2 diffusion 
along  the  coating/metal  interface  will  control.  In  the  transition  region  0.5  <m  <2,  neither 
mechanism  will  dominate.  Figure  3 delineates  the  three  regions. 

Higher  values  of  a/R  result  in  higher  values  of  t]  for  the  same  modulus,  rj  is  larger 
because  for  the  same  R,  diffusion  distances  under  the  coating  are  shorter,  so  that  concentrations 
remain  closer  to  the  value  in  the  inner  solution.  This  is  shown  in  Figure  4 for  m = 1 where 
the  slope  of  the  dimensionless  concentration  profiles  is  less  at  higher  values  of  a/R.  In  Figure 
5 higher  values  of  m result  in  steeper  predicted  profiles.  The  second  curve  from  the  right  shows 
the  approximate  profile  predicted  by  the  model  for  the  experimental  conditions  employed. 

Since  0^5  and  a/R  are  known  Eqs.(16),  (20)  and  (21)  can  be  solved  iteratively  for  the 
values  of  D^.,  m and  rj,  yielding:  T)^=1A  x ICf*  cmVs,  m = 12.1,  and  r}=0.055.  The  outward 
diffusion  coefficient  (DJ  value  of  sodium  through  the  coating  from  the  interface  is  one  to  three 
orders  of  magnitude  higher  than  that  reported  in  the  literature  for  an  epoxy  coating  free  film 
(6,7).  One  of  the  possible  explanations  for  the  higher  diffusion  coefficient  of  Na"^  obtained  in 
this  study  is  due  the  alkaline  environment  produced  by  the  corrosion  process  at  the  coating/steel 
interface.  The  alkaline  environment  increased  the  internal  stresses  in  the  coatings,  which 
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induced  macroscopic  or  microscopic  fractures  in  the  coating  and  facilitated  the  transport  of  ions. 


m,  dimensionless  modulus 

Figure  3.  Plots  of  rj  versus  modulus  m.  (See  nomenclature  section  for  definitions  of  tj  and  m) 


The  alkalinity  is  a product  of  the  corrosion  activity  under  the  coating  near  the  defects. 
This  comes  about  because  the  bare  metal  at  a large  defect  or  a scribe  mark  of  a polymer  coated 
steel  undergoes  anodic  reactions  (corrosion)  and  the  area  underneath  the  coating  surrounding  the 
defect  undergoes  cathodic  reactions  (8-10).  Sodium  ions  migrate  along  the  coating/steel  interface 
from  the  defect  to  the  cathodic  sites  (1,8,9)  to  neutralize  the  hydroxyl  ions,  which  are  produced 
by  the  cathodic  reactions,  forming  NaOH.  Thus,  the  liquid  at  the  coating/steel  interface  near 
the  defects  is  highly  alkaline.  Indeed,  pH  values  close  to  14  have  been  measured  at  these  sites 
(11).  The  alkalinity  of  the  liquid  at  the  coating/metal  interface  can  cause  the  disbondment  of  a 
coating  from  a metal  surface.  This  phenomenon  is  commonly  known  as  cathodic  delamination. 
As  stated  in  the  experimental  section,  we  observed  corrosion  at  the  bare  metal  in  the  inner 
cylinder.  It  is  likely  that  the  cathodic  delamination  occurred  underneath  the  coating  surrounding 
the  inner  cylinder.  This  was  evidenced  by  the  complete  detachment  of  the  film  from  the 
substrate  and  the  bright  appearance  of  the  substrate  after  the  film  was  removed  at  the  conclusion 
of  the  experiment.  This  may  provide  a partial  explanation  for  the  relatively  high  d value  derived 
from  the  model. 
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Figure  4.  Plots  of  dimensionless  concentration  0 versus  aspect  ratio  r/R  for  m = L 


Figure  5.  Plots  of  dimensionless  concentration  0 versus  aspect  ratio  r/R  for  m— 10. 
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The  alkaline  environment  at  the  coating/steel  interface  may  induce  internal  stresses  in  the 
epoxy  coatings  by  two  possible  mechanisms.  The  first  mechanism  was  the  development  of  an 
osmotic  pressure  gradient  between  the  outside  water  environment  and  the  coating/steel  interface. 
As  indicated  earlier,  the  area  underneath  the  coating  near  the  defect  underwent  cathodic  reactions 
and  this  area  is  rich  with  highly  hygroscopic  NaOH  corrosion  products.  As  the  result,  water 
activity  of  the  external  solution  is  higher  than  that  at  the  coating/steel  interface,  creating  an 
osmotic  pressure  difference  between  the  outside  and  the  inside.  Consequently,  water  flows  into 
the  coating  creating  a pressure  that  counterbalance  the  water  in-flow.  Inside  the  coating,  the 
pressure  generated  by  the  swelling  process  produces  an  elastic  stretching  of  the  polymer. 
Swelling  pressures  in  the  range  of  several  hundred  kgcm*^  have  been  reported  for  polyamide- 
cured  epoxy  films  when  they  absorb  water  (7).  Since  the  osmotic  pressure  gradient  is  a function 
of  the  water  activity  of  the  outside  environment,  the  use  of  distilled  water  in  the  outer  cylinder 
in  the  present  study  should  facilitate  the  coating  to  take  up  maximal  water  and  the  film  should 
stretch  to  its  greatest  extent. 

Another  mechanism  whereby  the  alkali  environment  at  the  coating/steel  interface  may 
induce  internal  stresses  in  the  coatings  is  through  an  increase  of  the  fixed  charges  in  the  coating. 
Epoxy  coatings  possess  fixed  charges,  which  are  produced  during  the  curing  process.  The  sign 
and  concentrations  of  these  fixed  charges  alter  substantially  when  the  environment  in  contact 
with  the  coatings  change  from  neutral  and  acidic  to  basic  (7).  Since  osmotic  pressure  is  directly 
related  to  the  concentration  of  fixed  charge  in  the  coatings  (7)  the  alkaline  environment  at  the 
coating/steel  interface  is  expected  to  increase  the  swelling  pressure  in  the  coating.  More 
hydrophilic  polymers  tend  to  swell  more  because  the  stronger  polar  or  ionic- water  interactions 
in  these  materials  will  cause  a greater  free  energy  decrease. 

The  existence  of  internal  stresses  in  swelled  polymers  and  their  effect  on  transport  have 
been  documented  (12,13).  In  the  presence  of  swelling  pressure,  the  diffusion  of  penetrant  has 
been  found  to  be  a linear  function  of  time  (instead  of  the  square  root  of  time  as  with  classical 
diffusion).  In  some  polymers,  the  internal  stresses  generated  during  diffusion  were  so  large  that 
the  films  cracked. 

The  relatively  high  value  of  m and  corresponding  low  value  of  77  both  indicate  that 
diffusion  along  the  coating/metal  interface  is  the  rate  limiting  step  for  the  transport  of  sodium 
ions  during  the  final  stage.  Note,  however,  that  the  ratio  of  the  diffusivity  in  the  channel  to  that 
in  the  coating  (Dg/DJ  is  quite  high  (from  23  to  230,  depending  on  the  value  of  5).  Diffusion 
is  effectively  two-dimensional  when  most  of  the  ions  get  to  pass  completely  underneath  the 
coating  before  they  start  to  diffuse  through  it.  In  general,  this  would  occur  with  a low  value  of 
m.  For  the  conditions  employed  in  the  present  experiment,  the  high  value  of  the  geometric 
factor  R^/h5  makes  m relatively  high  so  that  diffusion  will  not  be  two-dimensional.  This  can 
also  be  seen  in  Figure  5 where,  under  the  experimental  conditions  employed,  the  concentration 
is  predicted  to  fail  towards  zero  only  a short  distance  from  the  inner  cylinder.  The  steep 
concentration  gradients  indicate  that  the  principal  diffusional  resistance  occurs  in  the  channel 
under  the  coating. 
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6.  Summary  and  Conclusions 

Mathematical  models  were  developed  and  solved  for  the  two  stage  transport  of  ions  by 
diffusion  beneath  and  through  organic  coatings.  During  the  first  stage  transport  occurred  by 
radial  diffusion  in  the  channel  between  the  coating  and  the  metal  with  rapid  adsorption  on  the 
coating  underside  surface.  During  the  second  stage,  diffusion  also  occurred  outward  through 
the  coating.  Model  parameters  were  obtained  from  an  experiment  conducted  in  a specially 
designed  diffusion  cell  in  which  the  depletion  of  sodium  ions  in  a central  reservoir  was 
measured.  Model  predictions  of  concentration  versus  time  agreed  well  with  experimental 
results.  Most  of  the  sodium  ions  were  removed  during  the  first  stage,  as  indicated  by  the 
relatively  high  value  of  the  adsorption  parameter  7.  The  first  stage  occurred  more  rapidly  than 
the  second,  as  indicated  by  the  ratio  of  time  constants  for  the  two  periods.  High  values  of  the 
diffusion  modulus  m indicated  that  during  the  second  stage,  cation  transport  was  controlled  by 
diffusion  in  the  channel  between  the  coating  and  the  metal  substrate  rather  than  diffusion  through 
the  coating,  even  though  the  diffusivity  of  ions  under  the  coating  was  much  higher  than  that  for 
ions  through  the  coating. 
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7.  Nomenclature 


a:  radius  of  the  inner  cylinder. 

C^:  concentration  of  sodium  ions  in  channel  between  the  coating  and  the  metal  substrate. 
C^:  sodium  ion  concentration  within  the  inner  cylinder. 

C^o’.  initial  concentration  of  sodium  ion. 

C’ao^  value  of  Q at  t = 0;. 

Caw^  concentration  of  sodium  ions  in  outer  cylinder. 

D^:  effective  diffusivity  of  sodium  ions  through  the  coating  from  the  interface. 

Dgi  effective  diffusivity  of  sodium  ions  under  the  coating, 
f:  dimensionless  function  defined  by  Eq.(21). 

H:  height  of  cylinders, 
h:  coating  thickness, 

m:  dimensionless  modulus, 


m = 

Nfi  radial  flux  of  sodium  ions 

R:  radius  of  the  outer  cylinder, 

r:  radial  coordinate, 

rji  radius  of  the  adsorption  interface. 

Ar:  width  of  thin  shell  in  the  channel  located  at  radius  r. 
t:  time. 

V:  volume  of  solution  in  inner  cylinder. 


hh 
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V^:  volume  of  water  in  the  outer  cylinder, 
y:  dimensionless  radius,  r/a. 

w^:  molar  flow  rate  of  sodium  ions, 
a:  reciprocal  of  the  time  constant  for  model  I, 


6 D C 71  D 6 

a = ^ ^ 

V'i 

a’:  reciprocal  of  the  system  time  constant  for  model  11, 

! 2 Tz  a h r\ 

“ Vh 


7:  dimensionless  adsorption  parameter  in  model  I, 


Y = 


FC 


5:  effective  height  of  channel. 

Ps  : surface  molar  density  of  adsorbed  layer. 
7):  ion  penetration  parameter. 


2 af  hP^h 

\ 


[20] 


0:  dimensionless  concentration,  (CA-CAw)/(Ca(t)-CAw)- 
Gji  time  for  end  of  first  period  and  start  of  second. 

Ti'.  time  constant  for  model  I. 

Tni  time  constant  for  model  n. 
dimensionless  radius,  r/R. 
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Appendix  I:  Experimental  Data 

1.  Specimen  Characteristics 


Radius  of  inner  cylinder: 

a = 

0.6  cm 

Radius  of  outer  cylinder: 

R = 

2.25  cm 

Thickness  of  coating: 

h = 

150  ^m  = 

Height  of  solutions  in  cylinders: 

H = 

3 cm 

Volume  of  NaCl  solution  in  inner  cylinder: 

V = 

3.39  ml 

Volume  of  distilled  water  in  outer  cylinder: 

Vw  = 

44.23  ml 

2.  Concentration  CJppm)  of  Sodium  Ion  within  Inner  Cylinder 


= 9,000  ppm 


Time  (h) 

Ca  (ppm) 

Time  (h) 

Ca  (ppm) 

Time  (h) 

Ca  (ppm) 

0 

9000 

29 

1650 

109 

863 

1 

8760 

34 

1520 

114 

830 

5 

7800 

39 

1400 

119 

814 

6 

7050 

44 

1310 

124 

804 

7 

6390 

49 

1210 

129 

774 

8 

5690 

54 

1170 

134 

744 

9 

5000 

59 

1160 

139 

730 

11 

4020 

64 

1150 

144 

718 

13 

3420 

69 

1120 

149 

721 

15 

2990 

74 

1090 

154 

724 

17 

2710 

79 

1020 

159 

718 

19 

2410 

84 

1000 

164 

724 

21 

2220 

89 

988 

169 

721 

23 

2000 

94 

954 

174 

721 

25 

1870 

99 

929 

27 

1730 

104 

880 

20 


Appendix  U Computer  Codes 

1.  Program  MODELI 

2.  Program  MODEL2 

3.  Program  LESTSQ 

4.  Program  FINDYI 

5.  Program  BACKPLOT 

6.  Program  ETAVSM 

7.  Program  FINDDC 


21 


oonooon  oooon  n ooooo  o o o n o 


PROGRAM  MODELl 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

C: 

C 

C 

C 

C 

C 

C 

c 

c 

c 

C: 

C 

C 


****************************************************************** 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


TITLE: 

AUTHOR : 

DATE  Coded: 

ORGANIZATION: 

LOCATION: 


* 

INT  vs . Time  Calculation  * 

* 

Joanne  Zhang  * 

* 

May  6,  1991  * 

* 

Chemical  Engineering  Department  * 
Bucknell  University  * 

Lewisburg,  PA  17837  * 

* 


****************************************************************** 


PURPOSE: 


This  program  calculates  the  integral  form  of  the  mathematical 
model  I by  guessing  the  model  constant  Geimma,  applying  the 
experimental  data.  The  Integral  data  as  a function  of  time  can  be 
obtained.  A subroutine  in  IMSL  software  package  is  employed  to 
compute  the  integral  form. 


GLOBAL  CONSTANTS: 


INTEGER 

N 

PARAMETER 

(N  = 46) 

REAL 

CO 

PARAMETER 

(CO  = 9000) 

REAL 

R,  RO 

PARAMETER 

(R  = 4.5,  RO  = 1 

GLOBAL  VARIABLES: 

REAL 

CA(N),  T(N) 

REAL 

Y(N),  INTY(N) 

REAL 

INTC(N),  FY(N) 

REAL 

GAMA 

LOCAL 

VARIABLES : 

INTEGER 

I RULE,  NOUT 

REAL 

A,  ABS,  B,  ERRABS 

REAL 

EXACT,  F,  RESULT 

INTRINSIC 

ABS 

EXTERNAL 

F,  QDAG,  UMACH 

INITIALIZATION: 

ERREST , 


ERROR,  ERRREL 


BODY: 


OPEN  (UNIT  = 17,  FILE  = 'model. in’) 
OPEN  (UNIT  = 18,  FILE  = 'int.2g076') 


GAMMAM  is  given  a guessing  value  between  0.01  to  0.076. 


GAMA  = 0.076 
PRINT*, 

PRINT*,  'GAMMA 
PRINT*,  ' T 


GAMA 

F(Y) 


Integrating  by  using  the  guessing  GAMMAM  and  IMSL 
package . 

DO  10  I = 1,  N 

READ(17,  *)  T(I),  CA(I) 

IF  (I  .LE.  21)  THEN 

Y(I)  = SQRT( (1-CA(I)/C0)/GAMA  + 1) 

A = 1.0 
B = Y(I) 

ERRABS  ==  0.0 
ERRREL  = 0.0001 

I RULE  = 2 

CALL  QDAG(F,  A,  B,  ERRABS,  ERRREL,  IRULE,  RESULT,  ERREST) 

INTY(I)  = RESULT 

A = 1.0 
B = R/RO 

ERRABS  =0.0 
ERRREL  = 0.0001 

IRULE  = 2 

CALL  QDAG(F,  A,  B,  ERRABS,  ERRREL,  IRULE,  RESITLT,  ERREST) 

INTC(I)  = RESULT 

FY(I)  = INTY(I)/INTC(I) 

Data  output . 

WRITE(18,  99)  T(I),  FY(I) 

WRITE(*,  99)  T(I),  FY(I) 

ELSE 
END  IF 

CONTINUE 

TERMINATION: 


FORMAT (IX,  F5.0,  2X,  F8.5) 


END 


REAL 


FUNCTION  F(X) 


REAL 

REAL 

INTRINSIC 


X,  GAMA 

EXP 

EXP 


GAMA  = 0.076 


F = X*L0G(X)/(1+GAMA-GAMA*X**2) 


RETURN 

END 


PROGRAM  MODEL2 
C 
C 

n ****************************************************************** 


c 

* 

* 

c 

* 

TITLE: 

Model  II  for  CA/CO  vs.  Time 

* 

c 

* 

* 

c 

* 

AUTHOR : 

Joanne  Zhang 

* 

c 

★ 

* 

c 

* 

DATE  Coded: 

May  5,  1991 

* 

c 

* 

* 

c 

* 

ORGANIZATION: 

Chemical  Engineering  Department 

* 

c 

* 

LOCATION: 

Bucknell  University 

•k 

c 

* 

Lewisburg,  PA  17837 

* 

c 

* 

★ 

c 

c 

**********************************************************!>:******* 

c 

C: 

PURPOSE: 

c 

c 

c 

This  program  calculates  the  ion  concentration  (dimensionless) 

c 

as  a function  of  time 

by  applying  the  mathematical  model 

II . 

c 

The  best  slope  of  the 

linearized  the  experimental  data  is 

obtained 

c 

from  this  program. 

c 

c 


C: 

CONSTANTS : 

C 

C 

INTEGER 

N 

PARAMETER 

(N  = 46) 

REAL 

COP 

PARAMETER 

(COP  = 1210) 

C 

C 

C: 

VARIABLES : 

c 

REAL 

CA(N),  T(N) 

REAL 

LN(N),  TT(N) 

REAL 

B,  SUMYT,  SUMTS 

c 

c 

C: 

INITIALIZATION: 

c 

c 


OPEN  (UNIT  = 16,  FILE  = 'model2 . out2 ' ) 
OPEN  (UNIT  = 17,  FILE  = 'model . in ' ) 

C 


PRINT*, 

PRINT*,  ' T LN(C0P/CA)  ' 

C 

C : Data  input . 

C 

DO  10  I = 1,  N 

READ  (17,  *)  T(I),  CA(I) 

10  CONTINUE 
C 

C:  Calculate  the  best  slope. 

C 

SUMYT  = 0 
SUMTS  = 0 


C 


ooooo  ono 


DO  20  J = 21,  N 
TT(J)  = T(J)  - 49 
LN(J)  = LOG(COP/CA(J) ) 
WRITE(16,  99)  TT{J),  LN(J) 
WRITE(*,  99)  TT(J),  LN(J) 

C 

SUMYT  = SUMYT+TT( J) *LN( J) 
SUMTS  = SUMTS+TT( J) **2 
20  CONTINUE 

B = SUMYT/SUMTS 

Data  output . 

PRINT*, 

PRINT*,  ' BEST  SLOPE=  ',  B 


TERMINATION: 


99  FORMAT  (IX,  F5.0,  2X,  F7 . 3 ) 
C 

END 


oon  ooo  ooo  o oonnono  nnooo  onoo 


PROGRAM  LESTSQ 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C: 

C 

C 

C 

C 

C 

C 

C 

C 

C: 

C 

C 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


TITLE: 

AUTHOR : 

DATE  Coded: 

ORGANIZATION: 

LOCATION: 


★ 


Find  Best  GAMA  by  S**2  * 

* 

Joanne  Zhang  * 

* 

May  16,  1991  * 

* 

Chemical  Engineering  Department  * 
Bucknell  University  * 

Lewisburg,  PA  17837  * 

* 


***★★★***★*****★**★****★******★***★**★*★*********★******★*****★**★ 


PURPOSE : 


The  calculation  of  this  program  is  to  obtain  the  best  value 
of  the  model  constant  (in  model  I)  Gamma  by  minimizing  the  sum 
of  the  deviatioof  the  theoretical  data  and  the  experimental 
data . 


CONSTANT : 


INTEGER  N 
PARAMETER  (N  = 8) 


VARIABLES : 


REAL  GAMA(N),  SLOPE(N) 

REAL  S(N) 


INITIALIZATION: 


BODY: 


OPEN  (UNIT  = 20,  FILE  = ' gslope . dat2 ' ) 

OPEN  (UNIT  = 11,  FILE  = 'gs.dat2') 

PRINT*, 

PRINT*,  ' GAMMA  S**2  ' 

Data  input . 

DO  10  I = 1,  N 

READ(20,*)  GAMA(I),  SLOPE(I) 

Calculate  sum  of  the  deviation,  s**2  for  each  slope. 
CALL  SSQ(I,  SLOPE(I),  S(I)) 

Data  output. 

WRITE(11,  99)  GAMA(I),  S(I) 

WRITE(*,  99)  GAMA(I),  S(I) 

10  CONTINUE 


C 


noon  o non 


TERMINATION: 

99  F0RMAT(1X,  F5.3,  2X,  F9.6) 
END 


SUBROUTINE 

SSQ(NN,  SL, 

REAL 

TP 

PARAMETER 

(TP  = 39) 

INTEGER 

NN 

REAL 

SL,  SS 

SS  = (1/TP  - SL)**2 

C 

RETURN 

END 


oo  oooooon  oonoo  o nnooo  ooooo 


PROGRAM  FINDYI 


****************************************************************** 


c 

★ 

* 

c 

* 

TITLE: 

Yi  vs.  Time  Calculation 

* 

c 

* 

* 

c 

* 

AUTHOR: 

Joanne  Zhang 

★ 

c 

* 

Hr 

c 

* 

DATE  Coded: 

May  21,  1991 

* 

c 

* 

* 

c 

* 

ORGANIZATION: 

Chemical  Engineering  Department 

Hr 

c 

* 

LOCATION: 

Bucknell  University 

Hr 

c 

* 

Lewisburg,  PA  17837 

* 

c 

* 

* 

c 

C 

C 

C: 

C 

C 

C 

C 

C 

C 

C 

C: 

C 

C 


PURPOSE : 


This  program  calculates  the  integral  form  of  in  model  I as 
a function  of  time  by  using  the  best  value  of  the  model  constant 
Gamma.  A subroutine  in  ISML  software  package  is  applied. 


GLOBAL  CONSTANTS: 


INTEGER 

N 

PARAMETER 

(N  = 19) 

REAL 

TP 

PARAMETER 

(TP  = 39) 

REAL 

R,  RO 

PARAMETER 

(R  = 4.5,  RO  = 1.: 

GLOBAL  VARIABLES: 

REAL 

CA(N) , T(N) 

REAL 

YI(N),  YGUESS(335 

REAL 

INTC,  RHS(N),  DIF 

REAL 

GAMA 

LOCAL  VARIABLES: 

INTEGER 

IRULE,  NOUT 

REAL 

A,  ABS,  B,  ERRABS 

REAL 

EXACT,  F,  RESULT 

INTRINSIC 

ABS 

EXTERNAL 

F,  QDAG,  UMACH 

INITIALIZATION: 

BODY: 


OPEN  (UNIT  = 17,  FILE  .=  'model. in') 
OPEN  (UNIT  = 28,  FILE  = 'yt.dat2') 


GAMMAM  is  given  a best  value  by  minimizing  s**2. 


uuuo  o u uu  uuu  u u u uu  ouuo  ouuuouuuuu 


c 


GAMA  =0.07 
PRINT*, 

PRINT*,  'GAMMA 
PRINT*,  ' T 


YI 


GAMA 


Integrating  by  using  the  best  GAMMAM  and  IMSL 
package . 

A = 1.0 
B = R/RO 

ERRABS  = 0.0 
ERRREL  = 0.0001 

IRULE  = 2 

CALL  QDAG(F,  A,  B,  ERRABS,  ERRREL,  IRULE,  RESULT,  ERREST) 
INTC  = RESULT 


DO  10  I = 1,  N 

READ(17,  *)  T(I),  CA(I) 

RHS(I)  = INTC*T(I)/TP 

Find  Yi  at  each  T. 

DO  20  J = 1,  330 

YGUESS(J)  = (1.2  + 0.01*(J-1) )/R0 

A = 1.0 
B = YGUESS(J) 

ERRABS  = 0.0 
ERRREL  = 0.0001 

IRULE  = 2 

CALL  QDAG(F,  A,  B,  ERRABS,  ERRREL,  IRULE,  RESULT,  ERREST) 

DIF  = RHS(I)  - RESULT 

IF  (ABS(DIF)  .LT.  0.20)  THEN 
YI(I)  = YGUESS(J) 

GO  TO  50 


ELSE 

ENDIF 


20 


CONTINUE 


Data  output . 


50  WRITE(28,  99)  T(I),  YI(I) 

WRITE(*,  99)  T(I),  YI(I) 


CONTINUE 


TERMINATION: 


99  FORMAT (IX,  F5.0,  2X,  F8.5) 


END 


REAL 


FUNCTION  F(X) 


C 


REAL 

REAL 


X,  GAMA 
EXP 


INTRINSIC  EXP 


C 


GAMA  =0.07 


mM>W->.  t} 

F = X*L0G(X)/(1+GAMA-GAMA*X**2) 
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PROGRAM  BACKPLOT 


C 

C 

p •kic’kicicicieic’k'ic’kicic'k'klcicicicicicicfc'tcic'k'kic’kicic'kicicieicicicieic'tcicie’kic'k'kicif'kic'kic’kicicicicicicic'kie'k'k'k 


c 

* 

* 

c 

* 

TITLE: 

Backplot  for  Data  Set 

* 

c 

* 

★ 

c 

* 

AUTHOR : 

Joanne  Zhang 

* 

c 

* 

★ 

c 

★ 

DATE  Coded: 

May  22,  1991 

* 

c 

★ 

* 

c 

* 

ORGANIZATION: 

Chemical  Engineering  Department 

fc 

c 

* 

LOCATION: 

Bucknell  University 

* 

c 

* 

Lewisburg,  PA  17837 

"k 

c 

* 

k 

****************************************************************** 


c 

c 

C:  PURPOSE: 

C 

C 

C This  program  calculates  the  theoretical  data  points  by  applying 

C model  I and  model  I . The  best  value  of  model  constant  Gamma  and 

C the  best  value  of  the  slope  are  used.  The  set  of  data  will  be 

C employed  to  make  a backplot  of  the  experimental  data. 

C 

C 

C:  CONSTANTS: 

C 

C 

INTEGER 
PARAMETER 
REAL 

PARAMETER 
REAL 

PARAMETER 
REAL 

PARAMETER 


VARIABLES : 


REAL  T(N),  CA(N) 

REAL  yi(19) 


INITIALIZATION: 


100 

C 


BODY: 


OPEN  (UNIT  = 36, 

FILE  = 

' yt . dat2 ' ) 

OPEN  (UNIT  = 37, 

FILE  = 

' model .in') 

OPEN  (UNIT  = 38, 

FILE  = 

' bp2 . dat ' ) 

PRINT*, 

PRINT*,  ' T 

CA  • 

Data  input . 

DO  100  1=1,  19 

READ(36,  *)  T(I),  YI(I) 
CONTINUE 


N 

(N  = 46) 

CO 

(CO  = 9000) 

TP,  CP 

(TP  = 39,  CP  = 1400) 

GAMA,  ALPHA 

(GAMA  = 0.07,  ALPHA  = 6.05E-03) 


DO  10  I = 1,  N 


onnoo  ooo  ooo 


READ(37,  *)  T(I) 


Calculate  CA  as 


a function  of  T by  appliying  Yi  data. 


IF  (I  .LT. 

CA(I)  = 
ELSE 

CA(I)  = 
ENDIF 


19)  THEN 

C0*(1  - GAMA*(YI(I)**2-1) ) 
CP*EXP(-ALPHA*(T(I)-TP) ) 


Data  output. 

WRITE(38,  99)  T(I),  CA(I) 
WRITE(*,  99)  T(I),  CA(I) 

10  CONTINUE 


TERMINATION: 


99  F0RMAT(1X,  F5.0,  2X,  F8.1) 
C 


END 


oononnnon  ooooo  n nonn  n o o n n o n n n o n o n o n n n o n n o n o o o o o o n n 


PROGRAM  ETAVSM 


* 

* 

TITLE: 

Eta  vs . M Calculation 

•k 

* 

* 

ic 

AUTHOR: 

Joanne  Zhang 

* 

k 

* 

ic 

DATE  Coded: 

July  15,  1991 

k 

k 

* 

ORGANIZATION: 

Chemical  Engineering  Department 

k 

•k 

LOCATION: 

Bucknell  University 

* 

* 

* 

Lewisburg,  PA  17837 

k 

k 

★★★★★★★★★★★★★★★★★★★★★★★It****************************************** 


PURPOSE : 


This  program  calculates  the  the  effectiveness  factor  Eta  and 
as  a function  of  parameter  m with  the  ratio  of  a/R  as  a 
pareimeter . Bessel's  function  subroutines  in  the  IMSL  software 


package 

are  uesd  in 

CONSTANTS : 

REAL 

AR 

PARAMETER 

(AR  = 

0.5) 

INTEGER 

N 

PARAMETER 

(N  = 

101) 

VARIABLES : 

REAL 

EE(N) 

, M(N) 

REAL 

MM 

REAL 

F(N) 

REAL 

BSIO, 

BSIl, 

EXTERNAL 

BSIO, 

BSIl, 

INITIALIZATION: 

BODY: 


OPEN  (UNIT  = 22,  FILE  = 'em.kS') 

C 

MM  = 0.1 

DO  10  I = 1,  N 

IF  (I  .EQ.  1)  THEN 
M(I)  = 0.001 
ELSE 

M(I)  = (I-1)*MM 
ENDIF 
C 

F(I)  = (BSI1(M(I) ) *BSK1(M(I) *AR)-BSK1(M(I) ) *BSI1(M(I) *AR) ) 
+ /(BSI1(M(I) )*BSK0(M(I)*AR)+BSK1{M(I) ) *BSI0 (M( I ) *AR) ) 


oonoo  o non 


C 


EE(I)  = 2*F{I) *AR/(M(I) *(1-AR**2) ) 
10  CONTINUE 

: Data  output . 

PRINT*, 

PRINT*,  ' M EE  ' 

DO  30  I = 1,  N 

V7RITE(*,  55)  M(I),  EE(I) 

WRITE(22,  55)  M(I),  EE(I) 

30  CONTINUE 


TERMINATION: 


55  FORMAT (2X,  F7.3,  2X,  F7 . 5 ) 
C 


END 


ooonnono  onoon  n ooooo  o onoooooononoooooonononnoonnnno 


PROGRAM  FINDC 


* * 


* 

TITLE: 

M,  Eta,  DC  Calculation 

* 

if 

* 

if 

AUTHOR : 

Joanne  Zhang 

* 

if 

* 

if 

DATE  Coded: 

July  15,  1991 

* 

if 

* 

ORGANIZATION: 

Chemical  Engineering  Department 

* 

lit 

LOCATION: 

Bucknell  University 

* 

* 

ic 

Lewisburg,  PA  17837 

* 

★ 

**1t*********************************<r*******ifc*************1Ht****** 


PURPOSE: 


This  prograjn  calculates  the  prarameter  m in  f function  by 
guessing  m and  solving  equations ( 21) , (22),  and  (30)  simul- 
taneously to  find  the  best  value  of  m.  The  values  of  f and 
Dc  can  also  be  obtained.  Four  special  function  in  the  IMSL 
software  package  are  employed  in  this  program. 


CONSTANTS : 

REAL 

GAMMA 

PARAMETER 

(GAMMA  = 0.06) 

REAL 

ALPHA 

PARAMETER 

(ALPHA  = 9.5E-03) 

REAL 

ALPHAP 

PARAMETER 

(ALPHAP  = 5.0E-03) 

REAL 

AR 

PARAMETER 

(AR  = 0.267) 

REAL 

A,  R 

PARAMETER 

(A  = 0.6,  R = 2.25) 

REAL 

V,  H 

PARAMETER 

(V  = 3.39,  H = 1.5E-02 

REAL 

PI 

PARAMETER 

(PI  = 3.14) 

VARIABLES : 

REAL 

M,  ETA,  DC 

REAL 

DEDEL,  CM,  F 

REAL 

MG,  FG,  DCG 

REAL 

TOL 

REAL 

BSIO,  BSIl,  BSKO, 

BSKl 

EXTERNAL 

BSIO,  BSIl,  BSKO, 

BSKl 

INITIALIZATION: 

BODY: 


Calculate  DEDEL  by  using  ALPHA,  V,  and  GAMMA. 


ooono  ooo  o oono  o o non  n nnonnn  nnnn 


C 


DEDEL  = ALPHA/3 60 0*V*GAMMA/PI 


PRINT*, 

PRINT*,  ' DEDEL  = ',  DEDEL 

Calculate  CM  by  using  R,  DEDEL, and  H. 

CM  = R/SQRT( DEDEL *H) 

PRINT*, 

PRINT*,  ' CM  = ' , CM 

Get  frist  guessing  value  of  DC. 

DCG  = l.OE-11 

Calculate  M,  F,  then  check  DCG. 

DO  10  I = 1,  10000 
MG  = CM*SQRT(DCG) 

FG  = (BSI1(MG)*BSK1(MG*AR)-BSK1(MG)*BSI1(MG*AR) )/ 
+ (BSK1(MG)*BSI0(MG*AR)+BSI1(MG)*BSK0(MG*AR) ) 

DC  = (ALPHAP/3600*V/(2*PI*A*FG) ) **2*(H/DEDEL) 

TOL  = ABS(DC-DCG) 

IF  (TOL  .LE.  l.OE-12)  THEN 
GOTO  100 
ELSE 

DCG  = DCG  + l.OE-12 
END  IF 

10  CONTINUE 
100  CONTINUE 

PRINT*, 

PRINT*,  ' DC  = ' , DC 

Calculate  M and  ETA  by  using  DC. 

M = CM*SQRT(DC) 

F = {BSI1(M)*BSK1(M*AR)-BSK1(M)*BSI1(M*AR) )/ 

+ (BSKl(M) *BSI0(M*AR)+BSI1(M) *BSK0(M*AR) ) 

ETA  = 2*A*F/(R**2*(1-AR**2) ) * 

+ SQRT( (ALPHA/3600*GAMMA*V*H)/(PI*DC) ) 

Print  the  results  of  M and  ETA. 

PRINT*, 

PRINT*,  ' M = ' , M 
PRINT*,  ’ ETA  = ',  ETA 


TERMINATION: 


C 


END 
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